NBI-HE-00-10 



Macroscopic and Microscopic 
(Non-) Universality of Compact Support 
Random Matrix Theory 



G. AkemannW and G. Vernizzi^ 2 ) 

(1) Max- Planck- Institut fur Kernphysik 
Saupfercheckweg 1, D-69117 Heidelberg, Germany 

and 

(2) The Niels Bohr Institute 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 



Abstract 

A random matrix model with a a- model like constraint, the restricted 
trace ensemble (RTE), is solved in the large-n limit. In the macroscopic 
limit the smooth connected two-point resolvent G(z, w) is found to be non- 
universal, extending previous results from monomial to arbitrary polynomial 
potentials. Using loop equation techniques we give a closed though non- 
universal expression for G(z,w), which extends recursively to all higher k- 
point resolvents. These findings are in contrast to the usual unconstrained 
one-matrix model. However, in the microscopic large-n limit, which probes 
only correlations at distance of the mean level spacing, we are able to show 
that the constraint does not modify the universal sine-law. In the case of 
monomial potentials V(M) = M 2p , we provide a relation valid for finite-n 
between the fc-point correlation function of the RTE and the unconstrained 
model. In the microscopic large-n limit they coincide which proves the 
microscopic universality of RTEs. 



1 Introduction 



Random matrix models enjoy a wide range of applications in physics due 
to their property of being universal (for reviews see [jl], This property 
manifests itself in the independence of correlation functions from the choice 
of the distribution function V{M) ~ exp[— nTiV(M)], where M is an n x n 
matrix and V is a polynomial. Different classes of universality are found 
depending on the way the large-n limit is taken [||, [| and which part of 
the spectrum is investigated || Q, Q, §|- However, not in all applications 
a distribution function of the above form is realistic. Here, all the eigen- 
values (=energy levels) of the matrix M are coupled through the Jacobian 
after diagonalization. There are situations as for example in applications 
in Nuclear Physics, where in contrast to that the Hamiltonian of the model 
couples only few energy levels and is still very well described by the above 
random matrix model It is therefore very interesting to study deforma- 
tions and generalizations of the above distribution V(M) and to investigate 
first, if the correlations remain unchanged and second, if the property of uni- 
versality is maintained. In this work we study a deformation which preserves 
the symmetry of the matrix model, which will be the unitary transforma- 
tions of the Hermitian matrix M in our case. The symmetry of the model 
is directly related to the properties of the Hamiltonian under rotations and 
time-reversal There exists an interesting relation |J] between the re- 
stricted trace ensembles (RTEs) which we will consider and the so-called 
Wigner ensembles JlO| , pT| , where different matrix elements are weighted 
with different distribution functions, without being invariant under unitary 
transformations. 

It has been known only quite recently that there exist symmetry pre- 
serving deformations of the distribution function that destroy the prop- 
erty of universality |12|| . Two examples are the trace squared ensembles 
which were originally introduced in the context of Quantum Gravity [l4|] 
and the generalized RTEs which were introduced, in their simplest pure- 
quadratic form, by Bronk and Rosenzweig p5L [l6| in the context of Nuclear 
Physics. The non-universality does not necessarily spoil the applicability to 
physical systems. In fact the deformation may be introduced for physical 
reasons: indeed, trace-squared terms have been recognized as corresponding 
to higher order intrinsic curvature terms in the string action. Such terms 
have been added in order to cure an intrinsic instability of the theory related 



to a crumpled surface phase of the string world sheet |Tj|, 14 , 17 ] . Moreover 



in the framework of random surfaces, these additional terms are interpreted 
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as touching interaction terms which make the random surfaces touch each 
other. The generalized RTEs permit the same graphical interpretation as 
they have been shown to be a limiting case of a special trace squared ensem- 
ble []l8|]. The RTEs are defined such that the exponential weight function 
gets replaced by a constraint Vs(M) ~ 5(A 2 — 1/n TrV(M)) (or the 6 Heav- 
iside step function instead). By using the following representation of the 
5-function, 5(x) = lim^oo y/ir/l exp[— lx 2 ], the distribution Vs{M) can be 
easily seen to contain trace squared terms. The relation between RTEs and 
the trace squared ensembles was discussed in great detail for the spectral 
density in |[^] using saddle point techniques, and the canonical ensemble, 
the trace squared ensemble and the RTE were shown to agree to leading 
order. In also the two-point function was derived for the trace squared 
ensembles and shown to be non-universal using saddle point equations. Here, 
we will present a scheme to calculate general Appoint functions of the RTEs 
in different large-n regimes. 

Namely, one has to distinguish two different types of large-n limits - 
the macroscopic and microscopic limit - which may not all be affected by 
the deformation of the distribution function V(M). In the example of gen- 
eralized RTEs which we will study here we find that this is precisely the 
case. While the macroscopic universality is destroyed by the global con- 
straint, the microscopic correlations at short distances remain unchanged. 
In a sense, the canonical ensemble V(M) is replaced by its micro-canonical 
counterpart Vg(M). Therefore we have an explicit model of statistical me- 
chanics at hand, where the correlation functions of both ensembles can be 
calculated analytically and then be compared for discrepancies. Another 
remarkable property of RTEs is that they possess a finite support already 
at finite n. Similar to the canonical Gaussian matrix model the RTEs allow 
for an explicit calculation at finite re. This has been shown already in |l^] 
for the one-and two-point function and will be given here for general fe-point 
functions. Comparing the finite-n and re — > oo results the differences can 
be thought of as finite-size corrections, when interpreting n — > oo as the 
continuum limit. In the RTE these corrections appear in a different way 
than in the canonical model, due to the finite support at finite-n. Let us 
explain now in more detail how the large-n limit can be taken. 

In the macroscopic large-n limit no restrictions are made on the distance 
between different eigenvalues. This leads to smooth, universal two- and 
higher fc-point correlation functions for the canonical, unconstrained models 
||, 0]. Applications can be found in two-dimensional Quantum Gravity 
as well as the theory of transport properties of mesoscopic wires [Q] . For the 
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generalized RTEs we have shown in a previous publication together with 
our collaborators ||, that for a certain class of potentials V(M) the 2-point 
correlator is no longer universal. In this work we will extend these results 
to arbitrary polynomial potentials (see also |20[| ) and to all higher Appoint 
resolvents. Since in ref. Q it was also shown that to leading order all fc-point 
resolvents of the 5- and ^-measure are equivalent we will restrict ourselves 
here to the former one. 

In contrast to that, in the microscopic large-n limit correlations of eigen- 
values A, ft at a distance of the mean level spacing |A — fj,\ ~ 1/n are cal- 
culated. It is this kind of limit that finds a wide range of applications in 
nuclear physics, condensed matter physics (see e.g. Q) and has initiated ex- 
act analytical solutions in the study of Dirac spectra in QCD pi] . We will 
show that in generalized RTE with purely monomial potential V(M) = M 2p 
(which includes of course the quadratic case for p = 1) the constraint does 
not change the local properties at short distances, in the sense that the con- 
nected two-point correlator behaves according to the well-known "sine-law" 
of the Gaussian canonical ensemble ||. Indeed, this result holds also for 
higher fe-point correlation functions. This does not come as a surprise since 
a global constraint should not change the local properties. We believe that 
the same is true for more general potentials and that universality holds here 
as well. 

Consequently the present paper is split into two different parts. Section 
2 is devoted to the macroscopic large-n limit where we use loop equation 
techniques, closely following [22, 23]. Although the loop equations are origi- 
nally designed to calculate higher orders in the 1 / n-expansion we will restrict 
ourselves to the planar limit. Since we find non-universality to leading order 
in all connected correlation functions we do not calculate the likewise non- 
universal higher orders in 1/n. The difficulty to deal with the constraint will 
be treated in a similar way as the situation where the spectral density has a 
support consisting of several intervals p3|]. In these multi-band phases addi- 



tional constraints have to be imposed to make the solution unique [24 



We will restrict ourselves to hermitian matrices M only, for non-hermitian 
matrices see [20|. Since only the planar solution is needed for our results 
they can be easily extended to orthogonal and symplectic matrices using 
[p6| . Extensions to the complex matrix model are straightforward as well 
since the same loop equation techniques exist in the literature [|?j, . Very 
recently finite-n results have been obtained for the eigenvalue density of 
Gaussian ensembles with real symmetric and complex matrices psR . In the 
second part section 3 the microscopic large-n limit is investigated for the 



3 



RTE with a monomial potential V(M) = M 2p . Here, we generalize existing 
results for finite-n of previous publication [^] , by improving a technique al- 
ready used in ref . Q . Rescaling variables in the microscopic limit and using 
the inverse Laplace transform we are able to match the connected fc-point 
correlator to the well known "sine-law" behavior of the canonical ensemble, 
thereby proving the microscopic RTE universality. Let us stress that for the 
RTEs no orthogonal polynomial techniques are applicable. 

2 The macroscopic limit: non- universality 

In order to calculate all correlation functions for the constrained matrix 
model with an arbitrary potential V(M) = J27=i ^~M J we introduce an 
auxiliary potential W(M) inside the partition function^] 

Z s = JVM exp[-nTvW(M)] 5(n 2 A 2 - nTrV(M)) , (2.1) 

oo , 

W(M) = £ -^M j , (2.2) 
3=1 3 

where the two sets of variables {tj} and ({gj}, A) are taken to be inde- 
pendent. All Appoint resolvent operators can then be obtained by taking 
functional derivatives of Z$ with respect to W(p) as given below, where we 
then eventually set the auxiliary potential W to zero at the end. A similar 
trick has been used in ref. [12] in order to investigate a multi-trace random 
matrix ensemble fl3| , O], showing their non-universality as well. Further- 
more we use the complex representation of the <5-function to obtain 

Z s = J^JvMexp\ i -nTr(w(M) + ia(y(M)-A 2 )^. (2.3) 

If we had used instead the following representation of the ^-function, S(x) = 
lim^oo \/ir/l exp[— lx 2 ], we would have obtained the trace squared ensemble: 
Vi(M) ~ exp[2lnA 2 TiV(M) -ITrV(M) 2 ], with the strength of the touching 
interaction being proportional to I. However, it is not straightforward to 
derive and solve loop equations for such an ensemble. In particular, we 
cannot directly employ the non- universality results of |jl2|| , where it was 
crucial that the single- and multi-trace potentials were different. It is in the 



form eq. (2.3) that we can actually derive and solve the loop equations for 



: We added a trivial factor of n 2 inside the delta function. 
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the constrained model. Throughout the paper the same notation as in [^] 
is used, which is redisplayed here for completeness. The resolvent or 1-loop 
correlator is defined as 

OW = I( Tr ^-\ = i£<^ (2.4) 
yy; n\ p-M/g n^ Q p k+1 y ' 

and higher fc-point resolvents are given by 

G( Pl ,...,p k ) = n k - 2 (Tr^— . . . Tt-J— ) (2.5) 

\ pi-M Pk-M/ 5conn 

where conn stands for the connected part of the expectation value with 
respect to eq. (|2.3| ). They are defined such that the leading part is of the 
order 0(1). If we define the free energy Tg as follows 

Zg = exp[n 2 .F 5 ] (2.6) 

all resolvents can be obtained from it by successive applications of the loop 
insertion operator 

G( Pl ,..., Pk ) = ^(P^(P k -i)--~(Pi)^S + (2.8) 

In particular all higher resolvents can be derived from the 1-point resolvent 
alone. 

2.1 The loop equation in the planar limit 

The loop equation is derived in the usual way from the partition function 
eq. (^) by shifting variables M — > M + e/(p — M) and requiring it to be 
invariant under this shift, i.e. -^| e=0 = 0, 

/ *. KlM G(w) _ GW » + XM P )G [p ) , (2.9) 
jc 2-7TZ p — u n z aw 

where we have defined the effective potential 

V eff {M) = W(M) + iaV(M) . (2.10) 
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In eq.( |2.9D the integration contour C encircles the support [y,x] of the spec- 
tral density p(X) counterclockwise in the complex plane, not including the 
argument p £ [y,x]. The parameter ia inside the effective potential V e ff is 
determined by the constraint (TrV(M)) = nA 2 as a function of all coupling 
constants, as we will see in more detail below. This constraint can also 
obtained by requiring the invariance of Z$ under the shift a — > a + e. 

Let us stress again that the resolvents are given by differentiating with 
respect to W(p) and not the effective potential V e ff(p). Due to the a-integral 
in the partition function Zg eq. (|2.3| ) we have (TrM fc },5 ^ (iaTrM k )s which 
is needed to determine G(p) eq. ([OP- Furthermore let us note that eq. ( |2.9| ) 
looks almost identical to the loop equation of the unconstrained hermitian 
matrix model [^] defined in the next section eq. ( |3.1| ). However, the role 
of the above mentioned auxiliary potential W as well as the constraint will 
modify the results of f22]| . The constraint leads to similar complications as 
in the situation where the support of the spectral density consists of several 
intervals [p3fl . 

In order to solve the loop equation we introduce a 1/n 2 expansion for all 
fc-point resolvent operators (|2.5| ) 



G(pi,...,pk) = ^2-^-G g (pi,...,p k ) , (2.11) 

9=0 

where the leading part with genus g = (planar) is of the order 0(1). In ref. 
|J it was shown that for monomial potentials the expectation values of the 
RTEs possess such an expansion in 1/n 2 as well. Here we assume that the 
same holds true for all polynomials potentials. Inserting this expansion into 
the loop equation ( |2.9| ) and taking the large-n limit we obtain to leading 
order 

f ^Yk^l Go{uj) = Go{pf . (2 . 12) 

Jc 2iri p — to 

If we make the Ansatz that Go(p) has just one cut in the complex plane 
or equivalently the support of the eigenvalues consists of the single interval 
[y, x] we obtain 



Goip) = l(v: ff (p)-M(p)^(p-x)(p-y^j , 



(2.13) 



where the analytic function M(p) is given by 



M W - I , Sg£2- . (2.14) 

27rz (u) — p)y(ijj — x){uj — y) 
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For details of the derivation see for example ref. [23]. The final result can 
be written as follows, after deforming back the integration contour, 



1 I duj / (p-x)(p-y) 



G (p) = - & — e/JV - J-f — ^ . (2.15) 

2J c 2iri p-u V (u-x)(u-y) y ' 

To make the solution complete we still have to determine the endpoints x 
and y as well as the parameter ia in terms of the coupling constants of W 
and V and the parameter A. The first two equations can be obtained from 
the asymptotic behavior of Gq(p). According to the definition (|2.4D we have 



lim G(p) = - . (2.16) 

p— >oo p 

Since the leading term does not depend on n it comes from the planar part 
Gq(p) and we obtain the conditions 







1 I du V eff(u) 



2 Jc 2iri — x)(uj — y) 



1 = 1 / "IkM . (2.17) 

2 Jc 2iri (u — x)(uj — y) 

The third equation needed we obtain from the constraint on Tr^(M) which 
we rewrite in terms of the spectral density 

p(X) = — lim(G (A-ie) - G (X + ie 



2vr 

The constraint then reads 



2m e^o 

L M (\)y/(z-\)(\-v) . (2.18) 



A 2 = fd\ P {\)V{\) = f X ^M(X)yJ(x-X)(X-y) V(\) , (2.19) 



2tt 



which determines the parameter ia contained in Ai(X). This equation to- 
gether with the boundary conditions eq. ( |2.17| ) determines the planar resol- 



vent Gq(p) eq. (2.4) completely as a function of the coupling constants of 



W and V and of the parameter A. 
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2.2 Higher planar Appoint resolvents 

Starting from the the planar resolvent Gq{p) we can obtain all higher planar 
/c-point resolvents by successively applying the loop insertion operator 
to it, as it is given in eq. (|2.8| ). Here we use the fact that all resolvents 
have the same expansion in 1/n 2 eq. ( [HlD . For this purpose we introduce 



a set of new parameters M& and Jk, k £ N + . These moments usually play 
the role of universal parameters of the higher fc-point resolvents encoding 
all the information of the potential in addition to the endpoints x and y 
[p2|. We will then rewrite the loop insertion operator in terms of these new 



variables. This is done in order to make the successive application of 
to an algebraic procedure. Finally we calculate explicitly the non-universal 
planar 2-point resolvent Go(p,q) and comment on the general situation. 
Let us begin by defining 

Expanding the poles at x and y the moments can be explicitly written as 
functions of the coupling constants. Because of = ^ fc _ 1 - )! dX k-i M(A) 
and similarly for y the moments also characterize the multi-critical points 
of the model. We now rewrite the loop insertion operator eq. fl2.7|) in the 
following way 

d d . . dx d dy d dia d 

w <j>) = w (p) + w (r) Tx + w (r)^ + 

While the ^jy-(p) and (p) can be obtained in a straightforward way from 
the definition eq. ( J2.20 ) (see e.g. [f22j| ) the remaining unknown quantities 



are derived by applying (p) to the boundary conditions eqs. ( 2.17D and 



( |2.19| ) and solving a linear set of equations. This is done in the Appendix |A| 
with the result reading 



8 



j * w = _ ' (GpW -*,), 

alt p — y — y) 

' ''*%) = -4(G O (P)-0(P)) , (2.23) 



ia dW y ' B 



where 



B = ia f X dXV(X)— lim x (2.24) 
J y 2ni e->o 

G o (A-ie)-0(A-ie) - (G (A + ie) - <j>(\ + ie))' 

All quantities are expressed by elementary functions and the planar resolvent 
Go(p) eq. ( 2.15| ) and we have set already the auxiliary potential W = 0. 
The result for general W can be derived from eqs. (]Aj) and ([AH) . 

We are now ready to apply the loop insertion operator in the form 
( |2.21 ) to Gq(p) eq. (|2.15| ), which does not explicitly depend on the moments: 

G (p,q) = -^ 7 (q)G (p) 



dW 



2(p - q) 2 \4>(p) J A(q-p)(j)(p) 



^ dia . , 1 I duj V'(uj) j (p — x)(p — y) 



dW 2 Jc 2iri p — uj y (cj — — y) 

+ 4^)(^ Ml ^ (5) + ^ Jl ^ W )' ( ^ 5) 

after performing some contour integrals. If we set W = we can use the 
results eq. Q2.23] ) and finally obtain 

G (p,q) = Gl an {p,q) - L(G (p) - 4(p))(G Q (q) - tf>(q)) , (2.26) 

the planar connected 2-point resolvent of the constrained matrix model. The 
first part is the well known universal 2-point resolvent of the corresponding 
unconstrained or canonical matrix model (can) || 

{P,q) ~ <P-q) 2 \V(p-x)(p-y)(<i-x)(<i-y) ) [ ] 
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whereas the second part contains the non- universal terms Gq(p) and Go(q). 
Still, the result is given in closed form for an arbitrary polynomial potential 
V(M). Eq. ( |2.26j ) can be compared with the earlier result Q 



G (p, q) 



Gcan 




2p 



d p (pG (p))d q (qG (q)) 



(2.28) 



for the special case of monomial potentials V(M) = M 2p , p € N+. For p = 
1, 2 we have checked explicitly that the two results eq. ( 2.26 ) and eq. fl2.28| ) 
agree. The corresponding resolvents can be found in HQ and the parameter 
ia was already determined in [18|. From the procedure described above it is 
clear that also the higher fc-point resolvents will remain non-universal since 
the derivative ^{Pi)Go{Pj) always contains terms proportional to Go(pi) 
and G (j)j). 

Let us also briefly comment on higher genus contributions. Expanding 



the loop equation in 1/n together with eq. ( 2.11 ) one obtains for genus one 



(/C-2G (p))Gi(p) = G (p,p) 



where 



ICfip) 



(2.29) 



(2.30) 



fc 2ni p 

The right hand side of eq. (2.29) is easily obtained from eq. (2.26) by setting 
p = q. However, it is no longer a rational function in contrast to GQ an (p,p). 
Its non-universality will then translate to G\{p) after inverting the integral 
operator (IC — 2Go(p)) and thus to higher genera through 



9-1 

(it - 2G (p)) G g (p) = G g ,(p)G g ^ g ,(p) + -±-<j,)G g - X <p) , g > 1 . 

gl=l ' 



d 



dW 



(2-31) 

For this reason we do not go through the tedious procedure of finding a basis 
for (fC — 2Go(p)) now including also square roots and inverting it. 



2 In eq. (2.18) of ref. § the factor of l/2p is missing. 
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3 The microscopic limit: universality 



In this section we investigate correlations of eigenvalues at the distance of 
the mean level spacing D ~ 1/n, the so-called microscopic large-n limit. We 
will heavily exploit the knowledge about correlations in the unconstrained 
canonical model for finite as well as infinite n. Our main result is an ex- 
plicit relation between the canonical and RTE correlations for finite- ra which 
then serves to determine the microscopic RTE correlations and prove their 
universality. 

Let us recall the known results about the canonical ensemble which also 
fixes our notation. The partition function reads 

Z = J VM exp[-nTrV(M)] , (3.1) 

where V{M) is a polynomial. For the generalized RTE or micro-canonical 
ensemble 



Z s = J VM 5 (^A 2 - ■^TtV(M)^ , (3.2) 



where we do not need to introduce an auxiliary potential in contrast to the 
previous section. In order to obtain the same microscopic correlations for 
the two models we will eventually have to relate the coupling constants of 
the respective polynomial potentials V and V (as in the macroscopic limit, 
see [l^]). The Appoint density correlation function is defined as 

p(Ai, . . . , A fc ) = /-Tr 5(X 1 — M) ■ ■ ■ — Tr 5(\ k - M)\ , (3.3) 
\n n I 

and similarly for the delta-measure. Its connected part conn is related in 
the following way to the fc-point resolvent defined in the previous section eq. 

(PD 



p c (X 1 ,...,\k) = { -Tr *(Ai - M) ■ ■ ■ -Tr S(X k - M) ) . Ci.l) 
\ n n 

1 (±W 



n 



2k ~ 2 \27TiJ 



i™j ( n ai ) g ( Xi + aiie i • • • > A fc + 



a k ie) 



where we have not yet taken the large-n limit. In this section we deal with 
density correlations instead of resolvents because they can be given more 
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explicitly for finite-n. Namely in the canonical model all density correlators 
can be expressed in terms of the Kernel K n (X,p) of a set of orthonormal 
polynomials -P/(A) at finite-n p0[] 

p(Xi,...,X k ) = det [K n (Xi,Xj)] , (3.5) 

l<i,j<k 

p c {X x , . . . , A fc ) = (-l) fc+1 K n{M, X 2 )K n (X 2 , A 3 ) • • • K n (X k , A x ) , 

p 

where the sum is taken over the (k — 1)! distinct cyclic permutation P of 
the indices (1,2,..., k). The kernel and the polynomials are defined as 



i ii — ± 

K n (X,») = -e-m^yM)^P k (X)P k (p), (3.6) 

5 kl = [ dXe-^^PkimiX) . 



The use of the orthogonal polynomial method is only possible for the canoni- 
cal model because the measure exp[— nTrV(M)] inside the partition function 
eq. (|3.1| ) factorizes in terms of the eigenvalues Aj of the hermitian matrix 
M. For the RTE no such property holds which forces us to seek for other 
methods. Here we will make use of homogeneity properties for monomial 
potentials. 

Let us finally give the universal results for the canonical ensemble in the 
microscopic large-n limit, where we restrict ourselves to the origin scaling 
limit due to the local translational invariance of the canonical ensemble. As 
mentioned in the beginning we measure eigenvalues in units of the mean 
level spacing which is D = l/(np(0)) at the origin. Here p(0) is the mean 
eigenvalue density at zero, taken in the macroscopic large-n limit as given 
m eq. ( gig ). We then define new variables Z{ = Xi/D which are kept fixed 



in the large-n limit. Since D — > as n — > oo the variables Aj have to go to 
zero as well. In this particular limit the microscopic correlators are defined 

ps{z\,...,z k ) = lim (Dn) k p(ziD, . . . ,z k D) . (3.7) 

n — >oo 

This limit, which is well behaved and finite, can be investigated by using 
the Darboux-Christoffel formula for the kernel eq. (|3.6| ) and the asymptotic 



3 The factor n k appears as we had already defined the fc-point correlator in eq. (3.3) 
to be normalized to unity. The appropriate unfolding procedure is usually defined for 
un-normalized correlators (see e.g. which provides us with the correct pre- factor 

l/p(0) fc . 
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large-n behavior of the polynomials -Pfc(A) to obtain ||] 



lim DnK n (ziD, Z2D) 



sin(7r(zi - z 2 )) 



(3.8 



This is the universal sine-law which is valid for all polynomial potentials 
V(X). Together with eqs. ( ^ ) and (|3.7|) it completely determines all k- 
point density correlators, connected and not-connected, where at coinciding 
arguments we have lim^^oo DnK n (zD, zD) = 1, 



ps(zx,...,z k ) 



det 

l<i,j<k 



sin(7r(zj - Zj)) 
n(zi - zj) 



(3.9) 



and similarly for Pg(z\, . . . , Zkj- 

We conclude with the following remark. If we had taken the macroscopic 
large-n limit instead, the connected correlators p c (z\, . . . , z^) had been of the 
order 0(l/n 2fc ~ 2 ) as one can see from eq. ( [3.4| ) together with the fact that 
the (connected) resolvents are of order 0(1). However, when taking the 
microscopic limit keeping the Zi fixed, the asymptotic kernel eq. ( |3.8| ) is 
of order 0(1) and hence are the connected and not-connected microscopic 
Appoint correlators from eq. (3.5). Consequently, the knowledge of both 



connected and not-connected correlators, is equivalent here since they can 
be obtained from each other by adding or subtracting /-point correlators 
(I < k) of order 0(1). In that sense the microscopic large-n limit modifies 
the usual large-n factorization of correlation functions. 



3.1 Microscopic Appoint correlation functions 

As it has been mentioned already the correlation functions in the RTE 
can not be calculated using orthogonal polynomials because of the Dirac 8- 
function in the measure. However, in the particular case of purely monomial 
potentials V(M) = M 2p , the evaluation of the connected /c-point correlator 
in the microscopic limit is straightforward. We shall exploit some homogene- 
ity properties of this case, which make it possible to relate ensemble averages 
in the monomial RTE to ensemble averages of the same quantities in the 
canonical ensemble^] with potential V(M) = gM 2p . Let us define a function 
Fk{M; A) of the matrix M and of a set of parameters A = (Ai, . . . , A^), that 

4 This technique is a slight generalization of ref. jsjj, where only the macroscopic limit 
was investigated. 
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satisfies the following homogeneity property under a rescaling of the matrix 



F k (tM;X) = t a F k (M;t b X) for t, a, b G R . 



(3.10) 



A simple example for such a function is the operator Tr<5(A — M)/n which 
has a = b = — 1. Also any such product is a homogeneous function, as inside 
the average of the Appoint correlation function eq. fl3.3| ), having a = —k and 
b = — 1. In appendix |B] we derive the following formula for such homogeneous 
functions F k , which relates their canonical and RTE average 



(F k (M;X)) s 



Gf» a )*r(£) 



AT' 2 



(F k (M; ( 



t 



A)) 



t 2 p 



(A 2 ). (3.11) 



Here £ _1 [/i(t)](x) is the inverse Laplace transform of a function h(t), eval- 
uated at the point x > (for an integral representation see eq. (|B.5| )). Eq. 
( p. 11] ) holds for any finite n. If we choose the 1- or 2-point correlator from 
eq. ( |3.3| ) as an example we reproduce the finite-n results of ref. [18| for ps(X) 
and ps{\,jj) which were given in the case of a Gaussian potential p = 1. If 
we choose for general k to take (F k (M; X)}$ = ps(Xi, • • • , X k ), which obvi- 



ously fulfills the criterion ( 3.10 ), we obtain from eq. ( |3.11| ) the following 
expression for the /c-point RTE correlator 



ps(Xi, ■ ■ ■ , Afc) 



(W)^r(g) 

9 *-* 



A p 



-2 



p((a^)*Ai,...,(a^)*A* 



t 2 p 



(A 2 ) 



(3.12) 

Now we make use of the fact that at finite n the correlations p(Xi, ■ ■ ■ , X k ) 
from eq. (^^) can be written as a polynomial in all variables Aj times an 
exponential measure factor: 



p(Ai,...,A fc ) 



>f £ >) ,.x>>...xi £ 



2n-2 



(3.13) 



!i,...,ifc=0 

Due to this fact we can actually perform the inverse Laplace transformation 
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As(Ai, • • • , Afc) 



r — 

1 1 2p 



2n-2 c («) 



(P 2 ) 2 pAp 2 Zi,-,Zfc=0 (ffn 2 ) 2 ? 



x £ 



-i 



2j 



1=1 



fl ( A 2 - 1 , \ 2p \ 2n-2 >) \'i . . . \ l 

a [f 1 n l^i=\\ ) ^ C {h,...,l k } A l A h 



E 



(gn 2 ) 2 pA n p 2 h,...,i k =o {gn 2 ) 2 p 



■p / n 2 -k-Y, i l i 
1 V 2p 



1 



E A 



2p 



2p 



n f . 



, (3.14) 



where we have used the shift property £~ 1 [h(t)e at ](x) = C [h(t)](x + a), 
the linearity of the inverse Laplace transform and eq. ( [B.7D of Appendix 
B. Eq. ( |3. 14 ) is our first result, the finite-n fc-point correlation function for 
RTEs with monomial potential V{M) = M 2p in terms of the corresponding 
canonical correlator at finite-n with potential V(M) = gM 2p . In general, 
eqs. Q3.13 ) and ( 3.14| ) are different from each other at finite-n. This remains 
true in the macroscopic large-n limit, as we have seen in the previous section. 
In the remaining part we will show that in the microscopic large-n limit, 
however, they happen to coincide. 

In a first step, we determine the mean level spacing D$ = l/(np$(0)) in 
order to define the appropriate microscopic scaling limit. From eq. ( |3.14| ) 
at A = one can read off ps(0), since then the sum over all U collapses. We 
obtain 



ps(0) 



» 

'{0} 



(3.15) 



(gn 2 A 2 )^T^) (2pgA 2 )^ 

for its large-n value, with ojo} = p(0) being the macroscopic large-n limit 
of the canonical spectral density eq. ( |2.18| ) at the origin (which exists and 
is finite). In order to have the same mean level spacing as in the canonical 
ensemble we would have to set 2pgA 2 = 1. This identification of coupling 
constants occurs also in the macroscopic large-n limit in order to match the 
corresponding macroscopic spectral densities (see refs. |9[ [j^]). However, 
since we measure all correlations in units of D and D$ respectively we do not 
need to identify D = D$ since they drop out in the microscopic correlators 
anyway, as we will see below. 
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We will now take the microscopic limit analogue to eq. Q3.7| ) with rescal- 
ing Aj = Z{D$ of our finite-n relation ( 3.14| ). The large-n limit of the different 
factors can be obtained as follows, starting with the ^-function term in eq. 



e a 



2 2 pgA x ^ 

n 2p+l 



2p\ 



2 1 i—X 



6{A 2 ) 



(3.16) 



The remaining ^4-dependent terms yield 




n —k — ^ - 



2p 



» 2p+1 &W°}7 



2p 



+ 



( 1 + <^T 



(3.17) 



which still has to be evaluated under the sum over Zj's. Here we have used in 
the second step that the first factor inside the parenthesis is of 0(n 2 ) since 



J2k < fc(2n-2) , 



(3.18) 



i=l 



is at most of O(n). Because of p > 1 the corrections are sub-leading. The 
factor containing T-functions we evaluate together with the explicit factors 
of n in eq. ( |3.14| ) 



T? _ (k + Sjlj) 

. 2p 2p 



(3.19) 



( 9 n 2 )^ (fc+s ^ } r(| 

= (2 ra )-i(™)^l-^(A ; + ^Z i )(2p + fc + ^/ l ) + ---) • 

It remains to be shown that the second term and thus higher terms in the 
expansion are sub-leading. A naive counting from eq. ( |3.18|) suggests that 
this might not be the case. In order to use the microscopic results for the 
canonical correlators eq. (|3.9|) we put together our results obtained so far 



16 



lim (D s n) k ps(z 1 Ds, ... , z k D s ) = 

n — >oo 

- i im ( 2 P9^ 2 )^ 2 y c (n) I z 1 {2pgA 2 )^ Y _ / z k (2pgA 2 )h 

™ (c { o}) fc h ^ k J [h '- M \ n m ) \ n Ho} 

x (2pgA 2 )-^ (k+ ^ k) U ~ ^(k + J2l i )(2p + k + J2k) + .. ) 

\ " i=l i=l / 

(k k \ 

1 -4^(E^)(2p+E^) + -J 



2n-2 



x h S:-/ l ' ; '"(4)j ' (3 ' 20) 

Here, the pre-factors from eqs. ( |3.17| ) and ( |3.19|) have canceled with the 
factors of (2pgA 2 ) 2 p from the unfolding. Now we know from the canonical 
ensemble eqs. ( [3l3|) , Q and dU) that the limit of the 

sum over the L 

exists and is finite: 



1 2n ~ 2 ( z \ ll ( z \ lk 

ps{z\, ...,z k ) = lim -r V cf? j , — — • • • — — 

™ (c {0 }) fc kr ^ =0 l ' i, --'*> V nc io} / \ n m ) 

(3.21) 

Hence the term in eq. ( [3.20 ) proportional to l/(4pn 2 ) is indeed sub- leading 
and we have as a final result 

PS,s{zi 1 ...,z k ) = p s (zi,...,z k ) , (3.22) 



or more explicitly in terms of eq. (|3.9| ). Since our derivation holds for 
the RTE with an arbitrary monomial potential V(M) = M 2p , we have 
not only derived all k— point correlation functions but also proved their 
universality for the given class of potentials. Let us finally point out that the 



equivalence eq. ( |3,22| ) also holds for the corresponding connected A:— point 
correlation functions. As we have mentioned already at the end of the 
previous subsection, they are of the same order in the microscopic limit and 
they can be obtained from each other by adding or subtracting lower /-point 
correlators. 
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4 Conclusions 



We have shown for RTEs as an example of constrained random matrix mod- 
els that in one and the same model correlation functions may exhibit uni- 
versal and non-universal behavior in different large-n regimes. In particular, 
in the macroscopic large-n limit all planar connected /c-point resolvents are 
non- universal and a closed expression was given for the resolvent Gq(z,w) 
for an arbitrary potential. Hence when switching from the canonical to the 
RTE or micro-canonical ensemble, the delta-function constraint destroys the 
macroscopic universality. 

A different behavior appears in the study of correlations of eigenvalues 
at the scale of the mean level spacing 1/n. Here we recover the sine- law of 
the canonical ensemble and prove its universality for the class of RTEs with 
monomial potential. This leads us to conjecture that microscopic universal- 
ity holds also for more general RTEs. The result in the microscopic limit 
is not unexpected since a global constraint should not influence the local 
statistics of eigenvalues. 
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publication. Furthermore one of us (G.A.) wishes to thank the Physics De- 
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The work of G.V. is supported in part by MURST within the project of 
"Theoretical Physics of fundamental Interactions" . 



A The functions ^(p), ^(p) and §f (p) 

In this Appendix we apply the operator (p) eq. ( |2.21| ) to the boundary 
conditions eqs. (2.17) and (2.19) and solve the linear set of equations for 

ix /„\ dy / \ — j dia r„\ tt„;„„ j-t,„ ;j — j-:*.. 



the quantities g^(p), fy(p) and ffi (p). Using the identity 



d 
dW 



(A.l) 



we obtain 
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d P {p<P{p)) + \ 



„ , , . dy , . 

xMl ^ (p) + 2/Jl ^ W 



+ 



c 2ni 



from eq. ( 2,17] ) and 



C ( . , dx , . 

— Mi (p) 

2 V W w 



0(p) 
p — x 



(A.2) 



cLo V'(u 



p-yJ 

dw V'(u) 
c 2ni u — X 



1 

2 /c 2-7TZ p 




- y) 



(p-x)(p-y) 

from eq. ( 2,1S| ) after some calculation. Here we have introduced 

dX V(X) 



(A.3) 



C 



D 



2nx - X 
: dA F(A) 
2vry - A 



^/(x-AKA-y) , 
y/(x-\)(\-y). 



(A.4) 



We note that special care has to be taken when applying (p) to A4(X) 
which is then no longer an analytic function in A. Therefore the contour Coo 
has to be deformed in eq. ( p. 14 ) to contain only the pole and cut in the 
integrand, and not the new pole introduced by ^pp(p) (see also Appendix 
A in [ 23 1 ) . The linear set of equations eq. (A.2) and eq. ( |A,3 ) simplifies 
considerably when setting W = because of V e ff(uj) = iaV(u>) in that case. 
It is this limit which we will need in order give a closed final expression for 
the planar 2-point resolvent of the pure delta-measure without the auxiliary 
potential in eq. ( gjj| ). Using eq. ( gig) it follows 







„ , dx , . 



4>{p) 



p 



X 

2 

C 

2~ 



. , dx . , 
Mi (p) - 

Mi — (p) 



+ 
- x 

<P(p) 

p — X 

_ <p(p) 

p 
1 



+ 



.r 



0(g) 



<Kp) A 



+ 



p-yj 
<P(p)\ 
p-yJ 



2 ciia 
ia dW 



+^A 2 ^( p ) + ^(Go(p)-Hp)) 
ia dW ia 



(P) 



(A.5) 
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This can be easily solved for the desired quantities. We obtain 

j * w = m _ — — AGoip) - Hp)) 

aw p — y — y) 

1 r\ j 

:(P) = -o(Go(p)"^(p)) , (A.6) 



as given in eq. (2.23) where the following abbreviation has been introduced 



l2 . ~ C + D 



B = ia [ A + 2 . (A.7) 

V x -y ) 

One can easily convince oneself that it equals the form given in eq. ( |2.24| ). 

B RTE via inverse Laplace transform 

In this appendix we derive eq. ( |3.11| ) which expresses expectation val- 
ues with respect to the delta-measure in terms of averages with respect 
to the canonical measure eq. ( |3.1| ) using the inverse Laplace transform. Let 
F k (M; A) be a function of the Hermitian n x n matrix M and of the set of 
parameters A = (Ai, . . . , \ k ), such that it satisfies the homogeneity property 

F k (tM;\) = t a F k (M;t b X) (B.l) 

for some real t, a and b. In other words F k is a homogeneous function of 
degree a with respect to matrix elements My and of degree (—6) with respect 
to each parameter Aj. An example for such a function is the operator inside 
the average of eq. ( |3.3[) . The matrix integral considered in this appendix is: 



I[F k ] = J VMS f^A 2 - ^Tr[M 2 P]) F k (M; A) , (B.2) 

withp an integer number. It is proportional to the average (F k )s- Introduc- 
ing the complex representation of the delta function 2ir5(x) = fdyexp[ixy] 
and then scaling all matrix elements by a factor (gn 2 /(iy + + )) 1 / 2p , eq. 
O reads: 



„2 



(B.3) 
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where we have interchanged integrals. By using the homogeneity property 
( |B.1| ) we can rewrite the matrix integral as a canonical ensemble average, 



up to the normalization factor Z eq. (3.1): 



I[F k 



" d V JyA 2 



-oo 2-7T 



n +a 

Z (gn 2 )— CT X 



gn* 



2p 



iy + 0+ 



(F k (M; 



9n 



_b_ 



(F k (M; (*f Y* A)) 



t 2 p 



\iy + 0+ 



A)) 
(B.4) 



In the last step we used the complex representation of the inverse Laplace 
transform, i.e. 



1 ^+ioo+0 + 

£- L [h(t)]{x) = — dt e tx h{t) , x > 

2ni J-i 00+ o+ 



(B.5) 



In order to obtain the correct normalization of the delta-average we evaluate 



eq. ( B.4 ) in the case F = 1 (with a = b = 0): 



J[l] = (gn' 



l 2p 



(B.6) 



where we used the formula 
1 



c- 1 



(x) 



r( 7 + i 



-9(x) , Re(7) > -1 



(B.7) 



Finally, by putting together eqs. ( |B.4[) and flB.q) , the ensemble average 
(F k ) s = I[F k ]/I[l] with respect to the measure 5(A 2 - Tr[M 2 P]/n) reads: 



(gn 2 )2p ^ ( n 2 \ „ , 



A' 



2p 



(F fc (M;(2f) 2p A)) 



t 2 p 



(A 2 ) , (B. 



which is just eq. ( 3.11| ). 
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